NASA TM X-1394 


A QUASI-THREE-DIMENSIONAL METHOD FOR CALCULATING BLADE 

SURFACE VELOCITIES FOR AN AXIAL FLOW TURBINE BLADE 

By Theodore Katsanis and Lois T. Dellner 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 - CFSTI price $3.00 




A QUASI-THREE-DIMENSIONAL METHOD FOR CALCULATING BLADE 
SURFACE VELOCITIES FOR AN AXIAL FLOW TURBINE BLADE 
by Theodore Katsanis and Lois T. Dellner 

Lewis Research Center 
Cleveland, Ohio 

SUMMARY 

The aerodynamic design of turbine stator or rotor blades requires accurate deter- 
mination of the velocity distribution on the blades. For axial flow turbines, a quasi- 
three -dimensional compressible flow analysis has been used successfully for many years. 
The analysis is based on two methods - one to obtain a blade -to -blade variation of veloc- 
ity, and the other to obtain the radial variation in velocity. A combination of the methods 
gives a quasi-three -dimensional compressible flow analysis. Most of the calculations 
have been incorporated into a program. This report gives a description of the quasi- 
three -dimensional flow analysis and the FORTRAN IV computer program. 


INTRODUCTION 

The aerodynamic design of axial flow nozzles and turbine blades requires accurate 
determination of the velocity distribution on the blades. In the flow analysis, three- 
dimensional flow effects are of importance. Methods for accomplishing this type of anal- 
ysis for blade designs of medium or high solidity have been developed at Lewis Research 
Center (refs. 1 to 4). 

Reference 1 gives a method for determining the variation of velocity from hub to 
tip on a midchannel stream surface, and reference 2 gives a method for determining the 
blade -to -blade velocity variation. A combination of the two methods gives a quasi- 
three -dimensional flow analysis. Weight flow calculations are then based on a two- 
dimensional integration across a passage cross section. Most of the basic calculations 
needed for this analysis have been incorporated into a computer program to facilitate 
practical turbine design. The program is referred to as CTTD (Compressor and Turbine 
Division turbine design program). 


The purpose of this report is to give a complete description of the quasi-three- 
dimensional flow analysis and of the use of the CTTD computer program. This technique 
has been used at Lewis for 14 years, and the program has long been available to industry 
for their own use in turbine design. The original work of developing the analysis method 
and the CTTD program was done by the staff of the Compressor and Turbine Division at 
what was then the NACA Lewis Flight Propulsion Laboratory. 


SYMBOLS 


a parameter, eq. (Al) and (A2) 
b parameter, eq. (Al) and (A2) 

C curvature of streamlines on blade -to -blade surface, 1/ft 

2 

g gravitational constant, ft/sec 
m distance along meridional streamline, ft 
n distance along orthogonal to streamline, ft 

n q distance along orthogonal to streamline on blade -to -blade surface from suction to 
pressure surface, ft 

p 

p pressure, lb/ft 
R gas constant, (ft)(lb)/(lb)(°R) 
r radius, ft 

r c radius of curvature of streamline in radial -axial plane, ft 
T temperature, °R 
V absolute velocity, ft/sec 

W relative velocity, ft/sec 

w weight flow, lb/sec 

z axial coordinate, ft 

a angle of meridional streamline with the axial direction, deg 

/3 angle of streamline on blade to blade surface with the axial direction, deg (see 

fig. 3) 

y specific heat ratio 

p density, lb/ft^ 

2 



4o rotational speed, rad/sec 

Superscripts: 

absolute total condition 
” relative total condition 

Subscripts: 
calc calculated 

cr critical velocity 

e exit 

giv given 

h hub 

i inlet 

m mean 

mid midchannel 

p pressure 

s suction 

t tip 

6 tangential direction 


METHOD AND ASSUMPTIONS 

The objective of the analysis method is to calculate the quasi-three -dimensional 
velocity distribution satisfying continuity at a given channel orthogonal surface (see 
fig. 1). The weight flow may be specified in the calculation, or the calculation may con- 
tinue until the maximum (choking) weight flow for that channel orthogonal surface is 
determined. The velocity variation on the orthogonal surface is calculated from radial 
equilibrium (ref. 1) and stream filament theory (ref. 2). When this is done for several 
orthogonal surfaces, a velocity distribution over the blade surface is determined. The 
velocity distribution can thus be obtained for the guided channel formed by the portion of 
the passage where the blade -to -blade orthogonals extend from suction to pressure sur- 
face. The guided channel will not cover the entire suction surface. To obtain a velocity 
on the uncovered portion of suction surface some method must be used to estimate the 
location of the pertinent stagnation streamline. Since these velocities depend critically 
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on the location of the stagnation streamline, unreliable results may be obtained on the • 
uncovered portion of the blade. 

The basic simplifying assumptions used in deriving the equations used are 

(1) The flow relative to the blade is steady. 

(2) The fluid is a perfect gas. 

(3) The fluid is a non viscous gas. 

(4) The fluid velocity has no radial component. (The projections of the streamline 
on a radial-axial plane are straight lines parallel to the axis. ) 

(5) The midchannel line is a streamline, hereinafter referred to as the midchannel 
streamline. 

(6) The gas has a constant entropy from hub to tip at the midchannel streamline at a 
fixed axial coordinate. This assumption is used in calculating the midchannel velocity 
variation (ref. 1). 

The gas has a constant entropy blade -to -blade along an orthogonal to the stream- 
lines. This assumption is used in calculating the blade -to -blade velocity variation 
(ref. 2). 

(7) A line connecting the midpoint of the hub, mean, and tip orthogonals at a fixed 
axial location in the channel is a straight radial line. This assumption is used in calcu- 
lating the midchannel velocity variation (ref. 1). 

(8) The meridional streamline curvature varies linearly from hub to tip. (From 
assumption (4), meridional streamline curvature would be zero; however, the effect of 
wall curvature on radial equilibrium may be considered. ) 

(9) There is free vortex velocity distribution at the inlet to the blade; that is, (rV^ 
is a constant from hub to tip. This assumption is used in calculating the midchannel 
velocity variation (ref. 1). 

(10) The inlet absolute total temperature T^ is uniform. This assumption is used in 
calculating the midchannel velocity variation (ref. 1). 

(11) An additional assumption is necessary to calculate the velocity variation from 
blade -to -blade along an orthogonal. An option is provided for this assumption in the 
program. The usual assumption is that there is a linear variation of streamline curva - 
ture along an orthogonal. An alternate assumption is that there is a linear variation of 
radius of curvature along an orthogonal. There is no particular reason why one assump- 
tion is preferred over the other. For high solidity blading, it makes very little difference 
which assumption is chosen. 

Another possibility is provided in the program. In this case the blade surface 
velocities are calculated as above based on the assumption of linear variation of either 
curvature or radius of curvature. Then for the weight flow calculation, density and 
velocity along the orthogonal are computed by assuming a linear variation of static pres- 
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* sure along the orthogonal. This is not consistent with the assumption of linear variation 
of either curvature or radius of curvature and should be used with caution. 

* 

FLOW CHANNEL LAYOUT AND ENGINEERING DATA 

The initial steps in the design, those involving the development of the inlet and out- 
let velocity diagrams, are not described (see ref. 5 for further information). It is 
assumed that these velocity diagrams have been obtained, as well as the basic operating 
conditions of design weight flow, number of blades, gas to be used, operating speed, and 
inlet and outlet stagnation temperature and pressure. From this information an initial 
blade shape can be drawn (ref. 6). The following steps are followed to obtain a blade 
surface velocity distribution: 

A. Along a section midway between the hub and tip, a cylindrical development of the 
blade channel should be accurately laid out to scale several times actual size. Figure 2 
shows a typical blade channel - mean section. (This section need not be at constant 
radius, but the meridional streamline angle a should be small enough so that 

cos a « 1.) The midchannel streamline is then drawn midway between the suction and 
pressure surfaces, as shown in figure 2. This procedure is repeated for the hub and tip 
of the blade. The axial coordinate locations of the hub, mean, and tip sections relative 
to each other must be established. For this, an axial coordinate reference line is 
specified on each of the three blade sections. The relative angular location 0 of the 
hub, mean, and tip sections is not considered; the effect is usually negligible. 

B. Any number of axial locations can now be chosen on the three midchannel stream- 
lines. The distance along the midchannel streamline between each axial location should 
be measured for use in specifying loss distribution. At any given axial location a curve 
is drawn through the midchannel streamline from the blade -to -blade so as to be ortho- 
gonal to each blade surface and to the midchannel streamline, as illustrated in figure 3. 
This is done at hub, mean, and tip. The corresponding orthogonals at the three radial 
blade stations are positioned so that the intersections of the midchannel streamlines 

and the orthogonals are located at the same axial coordinate. In general, the intersec- 
tions of the orthogonals with the suction and pressure surfaces at the three radial sec- 
tions will not be at the same axial coordinates. These three orthogonals determine a 
section through the blade passage over which the velocity variation will be determined 
and across which weight flow will be calculated. The total length of the orthogonal ng 
is required as input for the program. Also, the blade surface curvatures C g and Cp 
must be measured at the ends of the orthogonals. These curvatures must be measured 
very carefully. The angle (3 is the angle that the midchannel streamline makes with the 
axis, and is considered positive in the direction of rotation; that is, p is positive if the 
tangential component of the velocity is in the direction of rotation. 
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If the inner and outer walls are not straight, the inner and outer wall curvature 
l/r c should be measured at each station. This curvature is considered positive if the 
wall is concave upward (fig. 4). This curvature is assumed to vary linearly between t 
hub and tip. The hub and tip radii complete the geometrical data required at each 
station. 

C . Operating conditions must be specified. These include the gas specific heat 
ratio y, operating speed co, and design weight flow per channel, w. Also required are 
W cr and p” at the hub, mean, and tip for each orthogonal. The inlet relative velocity 
W cr can be calculated using Vq . from the inlet velocity diagram together with the 
inlet absolute stagnation temperature T.': 


where 


W 


cr 


V 2yRg rj<M 

r + 1 


2u>r.V„ . - (wr)' 

rpf 1 _ rpt _ l t7 ? X 


2yRg 


(y - 1) 


( 1 ) 


( 2 ) 


From assumptions (9) and (10) (p. 4), T.’ and r.V fl . are independent of radius, and 

1 1 C/ j 1 

equations (1) and (2) can be used to determine W cr at any point in the passage as a 
function of radius alone . 

The inlet total stagnation pressure can be calculated from 



y/(y-l) 


(3) 


The exit relative total pressure p” can be found in the same manner utilizing the exit 
velocity diagrams. To allow for losses, the difference (p!’ - p”) is distributed along 
the length of the midchannel streamline. Usually a linear variation is used. Finally, 

p"=j£- (4) 

RT M 

D. The information determined in steps B and C is the information for items (l2) 
to (42) and (44) on the input data sheet (fig. 5). One sheet is required for each axial 
station. The calculation at any one station is independent of any other. It will be noted 
that the input sheet is designed so that only the numerical values which change from sheet 
to sheet need to be supplied. 
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The output from the program includes the blade surface velocities at each orthogonal 
and the corresponding weight flow. These velocity calculations are made for the initial 
( estimated value of W/W cr at the midchannel streamline for the mean blade section. It 
is possible for the program to determine the velocities corresponding to any desired 
weight flow (less than choking) or to determine the choking weight flow. If the velocity 
distribution corresponding to the design weight flow or some other operating weight flow 
is not satisfactory, the blade design should be altered, and steps A to C repeated to de- 
termine the new velocity distribution. 


DESCRIPTION OF PROGRAM INPUT AND OUTPUT FOR SAMPLE PROBLEM 

The input data sheet is shown in figure 5. The quantities filled in are for a sample 
problem. The output for this sample problem will be presented subsequently. 


CTTD 


General Instructions for Filling Out Input Data Sheets 

One page must be filled out for each z (axial) station. The basic calculation at each 


station is from a value of (W/W ) 


cr'mid, m 


(hereinafter called x) to the corresponding 


weight flow (WpQjg) and velocity distributions. 

(T) These two lines are for problem identification and will be printed as the heading of 
all z stations. Thirty -nine characters are allowed on each line, and they may be 
alphabetic, numeric, special symbols, or blanks. It is customary to fill out these 
two lines on the first page only and write OMIT on these two lines on all successive 
pages. However, if a new heading is desired, a page which is blank except for 
(IT) = KR4 = 9 must precede the page on which the new(T) occurs. 


be assumed in the calcula- 


© Jx-(a) “ u 

tion of the ve 

0 Jy . /sameas jxj u Unear variation to feme a p S r ^ ure ) is to be used in the calcu- 


linear variation in fcdSsTcurvature) is to 
velocity distributions across the orthogonals. 


lation of the weight flow. 
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(?) JZ = 1 if only the calculation for the one value of (44) = x is wanted. 

In this case, (l^) = w giv is not used. 

JZ = 2 if choke conditions are wanted. (Solve for the value of x which gives 

maximum value of w ) In this case w is not used. 

calc giv 

If there is any value of x for which (14) = w giy is equal to w ca j c , 
there are always two such values. 


JZ = 3 if the lesser of these two solutions is wanted (subsonic). 


JZ = 4 if the greater of these two solutions is wanted (supersonic). 

(?) (12) y ratio of specific heats 

(13) w wheel speed (rad/sec) 

(14) w giy weight flow PER BLADE, (lb/sec) 

First page. ^ 

Fill out @ - (14) and ( 5 ) = KR4 = 1. 


Any successive page. 

If values to be entered in 


those in (12) — 


© 

© 

@ 


/Q\/ differ from \ 

'O'yare identical toy 

of the immediately preceding page, KR4 = ( q ). If KR4 = 0, write 

OMIT in (12) — (14) . (?) = KR4 = 9 may be entered on an otherwise 

blank page. See explanation under (T) . 

( 1/r ] reciprocal of radius of curvature at hub (ft - *) 

. c /h 

^L/r c j reciprocal of radius of curvature at tip (ft - *) 

r^ radius at hub (ft) 

r^ radius at tip (ft) 

First page. 

Fill out @ @ and (?) = KR5 = 1. 

Any successive page. 

K values to be entered in @ - @ (ar^Wenti^al to) in © ** 

on the immediately preceding page, KR5 = (q\ If KR5 = 0, write 
OMIT in (l5) — (?8) . ' ' 
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. ® • ® c p 

, © , ® W cr 
. ® P"W C r 


length of the orthogonal between suction and 
pressure surfaces at hub, mean, tip (ft) 

curvature at intersection of orthogonal with 
suction surface at hub, mean, tip (ft - *) 

curvature at intersection of orthogonal with 
pressure surface at hub, mean, tip (ft - *) 

relative critical velocity at hub, mean, tip 
( ft/sec) 


both must 
be positive 
(nonzero) if 
JX= 2 


weight flow parameter at hub, mean, tip 
(lb/ft^ sec) 


First page. 

Fill in @ - 
If values for 


— (28) and (29) - (33) repeat those of (19) — (23), 
write OMIT in (24) — (33) and (lO) = KR6 = 1. If values differ, enter 


them and (10) = KR6 = 3. 


Any successive page. 

If all values to be entered in (l§) - (33) are identical to those in 

(19) — (33) of the immediately preceding page, write OMIT in (19) — (33) 

and KR6 = 0. If not, follow instructions given above for the first page. 

/3 relative (rotor) or absolute (stator) flow angle measured from the axis, 
(degrees), positive in the direction of rotation, midchannel only 

First page. 

If /S is constant from hub to tip, fill in (34) only; write OMIT in 
(35) - @ and (ll) = KR7 = 1. ~ 

If /3 is specified only at hub, mean and tip, fill in (34), (35), and (36); 
write OMIT in @ - @ and (ll) = KR7 = 3. 

If /3 variation from hub to tip is specified (nine values), fill in 
(34) - (42) and (U) = KR7 = 9. 

Any successive page. 

If all values in (34) — (42) are identical to those of the preceding page, 
write OMIT in (§ 3 ) — (42) and (ll) = KR7 = 0. If not, follow instructions 
given above for first page. 
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^43) An eight digit numerical code to identify this page uniquely is entered 

here and will be printed out preceding the output for this z station 
(can be used to code for engineer, case, and z -station number). 



( W cr)mid, 


m 


critical velocity ratio (referred to as x immediately below). 


Computation of w calc (and the velocity distributions) for this value of 
x will be performed first. If 

(T) = JZ = 1, calculation stops, and input data for next station is read. 

(T) = JZ = 2, x is increased automatically and calculations are continued 
until choking weight flow has been found. 

©=JZ = 3, x is modified until subsonic solution is found. (Point at 
which w ca j c = w gjy)- Guess as close as possible. 

©=JZ = 4, x is modified until supersonic solution is found. Guess low. 


Description of Computer Output 


An example of the output obtained is given in figure 6. The first output is a listing 
of the information on the input data sheet. In the sample output these are identified with 
the numbers corresponding to the numbers on the input data sheet. If nine values of /3 
(equally spaced from hub to tip) are not supplied, the missing values are calculated and 
printed with the input values. The remaining output is the calculated quantities. The 
first line is W g /W m ^ and Wp/W^^, each at hub, mean, and tip. Next are nine 
values (from hub to tip as 1 < IK ^ 9) of the parameters a and b of equation (A2) , 
which are labelled LITTLE A and LITTLE B. ALPHA is equal to a(r t - r h )/16. M(IK) 

is the numerical approximation to f T a(£)d£ and N(IK) is b • PH, PM, and 

J r 

m /*£ 

-1 a(£)d£ 

/ r m 

b(£)e di; at the hub, mean, and 

*m 


r a(£)d£ 
m 

tip, respectively. QH, QM, and QT are the numerical approximations to e 

The relative velocities ratioed to the relative critical velocity at the hub, mean, and tip 

for suction surface, midchannel, and pressure surface are followed by the values of 
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A(I,K) (pW/p”W cr ). The first line is at the hub at eight equal intervals from the suction 

surface to the pressure surface. The second and third line of A(EK) are at the mean and 

tip, respectively. On the next line N is the number of the iteration for this orthogonal, 

X is the value of (W/W ) for this iteration, WT FLOW CALC is the correspond- 

ing, m 

ing calculated weight flow w calc , and D is the measure of error 1 - w ca i c / w gi v - K 
more than one iteration is done, the values down to PH, PM, PT, etc. , do not change; 
hence, for the second and following iterations only the lines following this are printed. 

A maximum of five iterations is performed. If the convergence criterion is not satisfied, 
the statement NLIMIT (5) HAS BEEN REACHED will be printed. 

The message CALCULATION OF A(IK) AT ST. NO. 152+5 IN CTTD MP ASKS FOR 
LOG OF NEG. REMAINING CALCULATIONS THIS ITERATION THEREFORE INVALID, 
is usually caused by an error in the input cards. If this is not the case, the blade design 
must be changed. 
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CTTD COMPUTER PROGRAM 


The program consists of the main program and the two subroutines PABC and 

VSUBX. Subroutine PABC calculates the coefficients A, B, and C of the parabola 

o 

y = Ax + Bx + C passing through three given points. Subroutine VSUBX uses equa- 
tion (A4) or (A5) to calculate the ratio of the velocity at any point on an orthogonal to 
the midchannel velocity. 


C C ITD MAIN PROGRAM 

DIMENSION ENO I3),CSI3),CPI3J .WCRI3) ,PwCR(3) * VS ( 3 ) . V P 1 3 ) * W M ( 3 > , 

ISt 3),BFTA(9).AI FI 9I.EMI9) .CAPNI9) ,A(9) ,B<9) • WS (3 > ,WP<3> 

001 FOUIVAL ENCEIWSI 1 ) .WSH ) ,<WS<?) .WSM) , ( WS I 3 ) , WS T ) , t WP I 1 ) , W PH I . ( W P ( 2 
1 ) • WPM ) . I WP ( 3 ) * WPT ) 

DATA NL IM, TL IM/5..CC10/ 

400 FORMAT (40H1 

401 FORMAT { 40H 

40? FORMAT I 711) 

404 FORMAT ( 3F9.5) 

406 FORMAT ( 4F9.5) 

408 FORMAT ( 5F9.5) 

410 FORMAT ( 1H0 . 2 I 4 . E 9 . 5 ) 

41? FORMAT ( 1 1 HO GAMMA = ,E12.5,11H OMFGA = .E12.5*12H WG1VEN = . 
1 F12.5) 

414 FORMAT ( 106H0 N ZERO C SUB S C SUB P 

1 W SUB CR RHO W CR 1 /R SUB C R) 

416 FORMAT I 1OH0 HLB .7E15.5) 

417 FORMAT ( IOH MFAN .5E15.5) 

418 FORMAT I IOH TIP .7E15.5) 

4?0 FORMAT I17H0 BETA) 

4?? FORMAT ( IHO, 17. I4.E15.5.3I4 •113,311///) 

430 F0RMAT1/6H0 V SH =, F 1 2 . 5 , 6H VSM = E12.5.6H VST= .fc 1 2 .5 ,6H VPH=,E12. 
15. 6H VPM=,F12.5,6H VP T = , E 1 2 . 5 , /45HG IK LITTLE A LITTLE 

2 B ALPHA /) 

432 FORMAT I I4.3F15.5) 

434 FORMAT (31H0 IK M(IK) NllK) /) 

436 FORMAT ( I4.2F1 5.5) 

438 FORMAT I 6HC P H = , F 1 2 . 5 . 6H PM=E12.5,6H PT=,E12.5,6H OH=,E12. 

15.6H QM = , F 12 • 5* 6H QT=,F12.5) 

440 FORMAT ( //9H0 WSH * •E12.5.9H WMH * ,E12.5,9H WPH = .E12.5./9 
1H WSM = •F12.5.9H WMM = .E12.5.9H WPM = ,E12.5,/9H WST = • 
?F 1 2 • 5t 9H WMT * .E12.5.9H WPT = ,E12.5//) 

44? FORMAT ( 8H A ( I K ) = • 9E 1 2 . 5 ) 

444 FORMAT < 7H0 N = »I1.16H X = WMM/WCR = ,E12.5,18H WT FLOW CALC 
1= • F 1 2 . 5 • 7H 0 = • F 1 2 • 5) 

006 READ! 5,400 
OOfi READ ( 5,401) 

010 READ t 5,40?)JX,JY.JZ,KR4,KR5.KR6,KR7 

012 IFIKR4) 25,14.25 

014 IFIKR5) 3C. 16.30 

016 IFIKR6) 18,20.18 

018 I F ( KR 6— 1 ) 42.34,4? 

020 IFIKR7) 22,66,2? 

022 I F ( K R 7— 1 ) 24.48,24 

024 I Ft KR 7—3 ) 64.56,64 

025 I F ( K R 4— 9)26.6,26 

026 READ ( 5 • 404 ) GAMMA • OMF GA • WG I V 
028 GO TO 14 

030 R FAD I 5.406)KCH.RCT.RH.RT 
032 GO TO 16 
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0 34 RFAD ( 5, 408 )ENC( n,CS(ll,CPIl) *WCR <1 ) • PWCR 111 
036 DO 38 K = 1 * 2 

FNO(K*l)=FNOIK) 

CS (K + n-CS (K) 

CP(K+1 ) =CP ( K ) 

Wf.R(K + l)=WCR(K I 
038 PWCR ( K ♦ 1 ) »P WC R ( K ) 

040 GO TO 20 
04/ ufj 44 1-1*3 

04 4 RFAD ( 5* 408)FNG( I ) *C S( I ) *C P ( I ) *WCR (I ) * P WCR ( I ) 

046 GO TO 20 

048 RFAD ( 5*404 )RETA ( l ) 

050 DO 52 IK = 1 * 8 

052 BfcT A( I K + 1 )=BETA ( IK) 

054 GO TO 66 

56 RFAD ( 5*404)RFTA( I) *BF TA ( 5 ) *BETA ( 9) 

58 DO 60 K = 1 * 5* 4 

Tf MP = ( RFTA(K+4)-BETA(K) ) /4.0 

59 DO 60 J* 1*3 
KJ=K+J 

BFTAIKJ ) =BE TA I K J- 1 ) +TEMP 

60 CONTINUF 
062 GO TO 66 

64 RFAD < 5* 404) (8FTA( IK) *IK=1 *9) 

066 READ ( 5* 410 ID1* ID2*X 
068 WRITE I 6 • 4C0 ) 

070 WRITE (6*401 ) 

072 WRITE (6*412) GAMMA • OMEGA • MG I V 
074 WRITE (6*414) 

076 WRITE (6*416)EN0(1)*CS(1)*CP(1) .WCR ( 1 ) * PWCR ( 1 ) . RCH* RH 
078 WRITE ( 6* 417)FN0( 2 ) *C S( 2) *CP(2 ) *WCR(2) *PWCR(2) 

080 WRITE (6.418)tN0(3)*CS(3),CP(3).WCR(3) *PWCR(3) ,RCT*RT 
082 WRITE (6*420) 

084 WRITE ( 6*416)BFTA( l ) .BE TA ( 2 ) *8E TA ( 3) 

086 WRITE ( 6*417)BETA (4)*BETA(5) *8ETA (6) 

88 WRITE ( 6, 41 8) RETA ( 7 ) *BE TA ( 8) *BETA ( 9) 

U90 WRITE 16*422 )ID1*ID2*X*JX*JY*JZ *KR4 *KR5 * KRfe * KR7 
100 PM = 0*0 


OM 

= 1 • 0 

MPSFT 

= 1 

NO 

* l 

KSTAR 

= 1 

N 

= 0 

XOL D 

= 0.0 

TWO0M 

= 2 • O+UMFGA 

Cl 

= GAMMA- l .0 

C? 

= GAMMA /Cl 

C 3 

= C 1 /( GAMMA +1 • 0 ) 

C4 

= 1.0/GAMMA 

C5 

= 1.0/Cl 

DR 

* RT-RH 

DR6 

* 0R/6.0 

DR 16 

= DR/16.0 

DR24 

* DR/24.0 

DRC. 

P*0.0 

= RCT-RCH 
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101 DO 104 1*1,3 

CALL VSUBX (ENOU).CS(I) .CPU ) .P.JX.VX) 

104 vs< i ) = vx 

106 P = 1.0 

108 DO 110 I = 1.3 

CALI VSUBX IFNIOU) ,CS( I ) .CPU ) .P.JX.VX) 

110 VP ( I ) = vx 

112 WRITE (6,43G)VS(1),VSI2),VS(3) » VP ( 1 ) ,VPt2) * VP ( 3 ) 

114 P = 0.0 
116 DO 120 IK = 1 . 9 
RK =RH+P*DR 

RCK =RC H+P *DRC 

SNB = S I N I BETA! IK )* .01745329) 

B( IK I =TWOOM*SNB 
SNR2 = SNB* SNB 
CSR2 =1.0- SNB 2 
AK =<CSB2*RCK)-( SNB2/RK) 

ALF( IK ) = AK* DR 1 6 

118 WRITE (6,432)IK.AK,B(IK),ALEtIK) 

120 P * P+.125 

12? FMI4) * -ALE I 5 )-ALF I 4 ) 

EM I 3 ) = EM ( 4 ) -ALE I 4 )-ALF I 3 ) 

EM ( 2 ) = EMI 3)-ALF(3)-ALF(2) 

FMIl) = EM C 2 )-ALF ( 2 )-ALE II) 

EMI 6 ) = +AI F I 5 ) +ALF I 6 ) 

EM ( 7 ) = EMI 6 ) +ALF ( 6 ) +ALF ( 7 ) 

FM ( 8 ) = EM ( 7 ) +A L F l 7 ) +A LF I 8 ) 

FM ( 9 ) = EMI 8 ) +ALF I 8 ) +ALF ( 9 ) 

FMI5) *0.0 

123 WRITE I 6,434) 

124 DO 126 IK =1,9 

CAPNI IK I =B ( IK ) /E XP I EM 1 1 K ) ) 

126 WRITE I 6,436) IK. EMI IK).CAPN(lK) 

128 PH=-0R24*ICAPNI 1) K. APN ( 5) +2. 0*C APN 1 3 ) *4. 0* I CAP M 2 ) +C APN 1 4) ) ) 
PT=+DR?4*<CAPNI5)+CAPN19)+2.0*CAPN( 7) +4 . 0* I CAP M 6 ) +CAPNI8) ) ) 
OH = FXP I FM I 1 ) ) 

OT =EXP I EM I 9 ) ) 

130 WRITE I 6, 43 8 )PH, PM . P T, OH • OM ,0T 
140 WM I 1 ) =QH*| X*WCR(2)-PH)/WCR(1) 

WMI3) =0 T*( X*WCR (?) — PT)/WCR(3) 

WM I 2 ) = X 

WSH =WM I 1 )* VS ( l ) 

W ST =WM( 3 ) *VSI 3 ) 

WPH =WM I 1)*VP<1 ) 

WPT =WM( 3)*VP< 3) 

WSM =X*VS<2) 

WPM =X*VPI2) 

142 WRITE I 6. 440 ) WSH, WM (1),WPH,WSM,WM(2) ,WPM,WST,WM(3) ,WPT 
144 DO 166 1=1, 3 
146 SIGA =0.0 
P=0.0 

ENOC =ENO< I) 

148 GO TO I 150. 150.158), JY 
150 CPC =CP( I) 

CSC =CSII) 

WMC =WM(II 



15? on 154 IK - 1 * 9 

CALL VSUBXC FNGC.CSC .CPC . P • JX . VX) 

ZK =VX*WMC 
TFMPO = l.-C3*ZK*ZK 
IFCTFMPOU 153*1153*1 53 
115 3 TFMPO = ABS( TEMPO I 
WR 1 f E i 6 « i 234 * 

12 34 FORMAT l 56H CALCULATION OF AUK) AT ST. NO. 152+5 IN CTTO MP ASKS 
12341. 60HF0R LOG OF NEG. REMAINING CALCULATIONS THIS ITERATION THERE t 
1234? 13HFORE INVALID. ) 

153 At IK ) = ( T EMPG**C 5 ) * ZK 
SIGA = S I GA +A t IK) 

154 P =P + . 125 
156 GO TO 164 

158 FLS =t 1.0-C3*WS( I>*WS<I) >**C2 
FLP = t 1.0-C3*WP< I )*WPU ) )**C2 
160 DO 1 62 IK = l • 9 

CAPJ =ELS+P*( ELP-EL S) 

TEMP =CAPJ**C4 

At IK) =TFMP*SQRT< ( 1 . 0-C A P J / TE MP ) /C3 ) 

SIGA = S IGA+A ( IK ) 

167 P =P+ • 125 

164 WRITE ( 6. 44? )( A UK) , IK = 1,9) 

TMP = . 5*t At l ) +A ( 9 ) ) 

166 S(l) =< SIGA-TMP)*(EN0C/8.0> 

168 N =N+1 

W = DR6*t St 1 )*PWCR(1) + 4.G*S(2) *PWC 8 1 2 ) + S ( 3 ) *PWCRI3) ) 

D = 1 • C- ( W/WGI V) 

1/0 WRITE (6.4441N, X.W.D 
17? GO TO C 10. 174. 174. 174) .JZ 
174 IF CNLIM-N) 176,176.180 
176 WRITF (6. 500)NLIM 

500 FORMAT I 13H0 NLIMIT C.I3.18H) HAS BEEN REACHED ) 

178 GO TO 10 

180 GO TO ( 182.186. 198, 198). JZ 
18? WRITE C6.501) 

501 FORMAT < 19H ERROR STOP AT 180 ) 

184 GO TO 10 

186 GO TO C 188.19?,256).KSTAR 

188 KSTAR = 2 

190 GO TO 710 

192 KSTAR =3 

194 GO TO ?40 

198 IF( ABS ( 0 )-~TL IM) 10.2CC.200 
200 GO TO ( 202. 230, 24?, 256). NC 
?0? NC = 2 

204 IF( A8SC D)— .05)206,210.210 
?06 OX =.05 
208 GO TO 212 
210 OX =.20 
212 X)=X 
W 1=W 
Dl = 0 

214 IF(D.GT.O.O) X=X+OX 
216 IF(D.LT.O.C) X=X~OX 
226 GO TO 140 
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230 NC= 3 

23? IFI 0*01)238. 234,240 
234 WRITF (6.503) 

50 3 FORMAT ( 19H FRROR STOP AT 232 ) 

236 GO TO 10 

238 DX =.5*0X 
MPSFT =2 

240 X? = X 
W2 = W 

02=n 

241 GO TO 214 

242 IF( 0*02)246, 244. 25G 
244 WRITF (6.504) 

504 FORMAT 1 19H FRROR STOP AT 242 ) 

246 GO TO 10 

248 MPSFT =2 

250 GO TO ( 252,256),MPSET 
252 X 1 = X 2 
W 1 = W 2 
01 = 02 
X 2 = X 
W2=W 
02 = 0 

254 GO TO 214 
256 X 3 = X 
W 3 = W 
03=0 

CALL PARC ( XI, X2 .X3.W1.W2. W3.APA8 * ft PAR .CPAB) 

260 GO TO ( 262.266.29C.290), JZ 
262 WRITE (6.505) 

505 FORMAT ( 19H FRROR STOP AT 260 ) 

264 GO TO 10 

266 X =-BP AB/ ( APAR + APAB ) 

268 IF( .001-ABS( X-XOLO ) )270.270.10 
270 XOI 0-X 

27? IF(Wl-W2) 274.280.280 
274 IF(W1-W3)276, 14C.14C 
276 X 1 =X3 
W 1 =W 3 
01 = 03 

278 GO TO 140 
280 I F( W2-W31282, 14C. 140 
282 X 2= X 3 
W 2 = W 3 
02=03 

284 GO TO 140 

290 DISC=RPAB**2-4. C*APAB* ( C PAB- WG IV) 

292 IF ( DISC ) 294*298*298 
294 WRITE ( *. 506) 

506 FORMAT ( 55H PARC FIT GIVES NEGATIVE DISCRIMINANT. PRUBA6LY CHOKED. ) 
296 GO TO 10 

298 X PL = ( - BP AB + SQR T ( 0 I SC ) ) / ( A PAB+ A PAB ) 

XM I =( -BPAB- SQR T ( 0 I SC ) ) / ( A PAB+ A PAB ) 

300 GO TO I 202.302.306.3081.37 
302 WRITF (6.507) 

507 FORMAT ( 19H FRROR STOP AT 300 ) 


304 GO TO 10 

506 X = AMINK XPL,XMl> 

GO TO 316 

308 X = AM AX l ( XPL « XM I ) 

316 IF <ABS(0i>-ABS<D2) >318, 320,320 
3IH IF i ABSI 02#- ABS 103 i ! 328 ? 3 22. 32 2 
320 IF (ABS(01)-ABS<03) >328.326.326 
322 X2 = X 3 
W2=W3 
02 = 03 

324 GO TO 328 
326 X 1 =X 3 
Wl = W3 
01=03 
328 NC=4 
330 GO TO 140 
END 


C SUBROUTINE PABC 

SUBROUTINE PABC < X 1 , X2 • X3 , Wl . W2 . W3 * A ,0 ,C ) 
C9 = X1 + X2 
C11=X1*X1 
C 1 5 = X 3- X 1 

Cl 8= <W2-Wl)/<X2-XU 

A = t C 1 5*C 1 8-W3+W1 )/(Ci 5*C 9— X3*X3+C 1 1 ) 

B=C18-C9*A 

C = W1-X1*8-C11*A 

R FT URN 

ENO 


C SUBROUTINE VSUBX 

SUBROUTINE VStJB X IEN*CStCP*P*JX*VXi 
700 GO TO 1 7C2 » 706 ) , J X 

702 VX=FXP( EN*(CS*< .5-P )♦. 6*<C S-CP) * < P*P-. 251) 1 
704 GO TO 714 

706 IF CCS-CP) 712*708,712 
700 V X = FXP I FN *C S* ( • 5- P 1 ) 

710 GO TO 714 

712 VX = < 2 . 0 *UP*(CS-CP>+CP)/<CSfCPm**<EN*CS*CP/<CP-CSn 
714 RETURN 
END 
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Program Variables in Main Program 


A 

array for nine values of pW/p"W cr 

AK 

a 

ALF 

array for nine values of a(r^. - r^)/16 

APAB 

O 

coefficient of x calculated by subroutine PABC 

B 

array for nine values of b 

BETA 

array, input, (3 

BPAB 

coefficient of x calculated by subroutine PABC 

C 1 

Y - 1 

C 2 

y/(y - 1) 

C 3 

(y - l)/(y + l) 

C 4 

l/y 

C 5 

i/(y - i) 

CAPJ 

p/p” 

CAPN 

array, b‘e -EM 

CP 

array, input, C 

Jr 

CPAB 

constant coefficient calculated by PABC 

CPC 

temporary storage for current value of CP 

CS 

array, input, C g 

CSB2 

cos 2 0 

CSC 

temporary storage for current value of CS 

D 

1 - w calc/ w giv 

D1,D2,D3 

temporary storage for various values of D 

DISC 

temporary storage 

DR 

r f r h 

DR16 

(r t - r h )/16 

DR24 

(r t - r h )/24 

DR6 

(r t - r h )/6 

DRC 

(1/ r c ) - (l/r c ) 
t h 
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DX 


arbitrary change in 
of weight flow 



ELP 

ELS 

EM 


ENO 

ENOC 

GAMMA 

i 

I 

I 

EDI 


ED2 


Pp/P" 

p s /p" 

arrav. r a d£ (trapezoidal rule is used) 

■ 

m 

array, input, ng 

temporary storage for current value of ENO 
input, y 

index of most loops executed three times (for hub, mean, tip) 

input, z station code, 4-digit number supplied as input and printed 
out for identification 

input, z station code, 4-digit number supplied as input and printed 
out for identification 


IK 

J 

JX 

JY 

JZ 

K 

KJ 

KR4, KR5, KR6, KR7 
KSTAR 


index of most loops executed nine times 

index of an input storage loop 

input 

input 

input 

index for two loops controlling storage of input 
temporary storage K + J 
input, read control switches 

switch - controls branching for successive values of x when 
choke solution is wanted; initialized by program and auto- 
matically stepped 


MPSET switch - controls branch to continue calculations or call subroutine 

PABC; initialized by program and altered as a result of calcu- 
lations 


N counter for number of times weight flow is calculated (When 

N = NLIM calculation for that page of input is stopped and a 
message written at end of output. ) 
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NC switch - controls branching for successive values of x when subsonic or 

supersonic solution is wanted 

NLIM constant (supplied in DATA statement) limiting number of times weight flow 
calculation can occur for 1 sheet of input data 

OMEGA input, w 

P n/n Q 


PH 

J h b • e" EM dr (Simpson’s 
r m 

PM 

0 

PT 

f * b • e~ E ^ dr (Simpson's 

Jr m 

PWCR 

array, input, p” W cr 

QH 

e EM at hub 

QM 

1.0 

QT 

EM , . . 
e at tip 

RCH 

input, (1/r ) 
c h 

RCK 

current value of l/r c 

RCT 

input, (1/r ) 

t 

RH 

input, r h 

RK 

current value of r 

RT 

input, r t 


array, values of 


r n o pw 
Jo p" W cr 


dn at hub, mean, and tip 


SIGA temporary storage, £A(K) 

SNB2 sin 2 /? 

SNB sin /? 

TEMP temporary storage 

TEMPO temporary storage 
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TLIM 


TMP 

TWOOM 

VP 

VS 

VX 

W 

W1,W2,W3 

WCH 

WGIV 

WM 

WMC 

WP 

WPH 

WPM 

WPT 

WS 

WSH 

WSM 

WST 

X 


X1,X2,X3 


constant (supplied in 
w 


TLIM > 


1 - 


calc 


w„ 


giv 

temporary storage 


DATA statement) (Calculation is terminated if 
whenever subsonic or supersonic solution is sought. 


) 


2o> 


array, three values of velocity (calculated in subroutine VSUBX) for 
pressure surface at hub, mean, and tip 

array, three values of velocity (calculated in subroutine VSUBX) for 
suction surface at hub, mean, and tip 


ratio of velocity to midchannel velocity 

weight flow calculated for current value of ( \ 

temporary storage for various values of weight flow 
array, input, W cr 


input, w giy 

w 

array, midchannel values of at hub, mean, and tip 

W cr 

current value of WM 


array, W/W cr on pressure surface at hub, mean, and tip 
WP(1) 

WP(2) 

WP(3) 

array, W/W cr on suction surface at hub, mean, and tip 
WS(1) 

WS(2) 

WS(3) 

input, also the current value of 


r_W 

m 


for which a weight flow is 


cr , 


mid, m 


calculated 

temporary storage for various values of x 
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XMI 


lesser of two values of x at which line y = w . intersects parabola 
0 

y = APAB&T) + BPAB(x) + CPAB 

XOLD temporary storage for a value of x when choking weight flow is sought 

XPL larger of two values of x at which line y = w^ iv intersects parabola 

y = APAB(x 2 ) + BPAB(x) + CPAB 

ZK current value of W/W at any point 

Li 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, March 17, 1967, 
720-03-01-35-33. 
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APPENDIX - ANALYTICAL EQUATIONS 


In reference 1 a differential equation is derived for the variation of velocities in a 
radial direction based on the assumption of a midchannel stream surface (axial symmetry) 
which is radial at any fixed axial position (i. e. , — = 0 on the midchannel stream sur- 
face): 9r 


dW 

an 


= aW - b 


where 


? 9 

_ eos*/3 _ sin /3 


r cos a 


b = sin /3 


2oj 

icos a 


■ + tan a 


dW ( 

dm 


(Al) 


The coordinate system is shown in figure 7. 

With the assumption that a is sufficiently small so that cos a is nearly one and 
tan a is nearly zero, we can set dn = dr to obtain 


where 


SW = aW - b 
dr 


2 2 
a = cos p _ sin /3 


r c r 


(A2) 


b = 2w sin j3 

An analytical solution of equation (A2) is 


L a(£)d£ 
W(r) = e m 





for the midchannel stream sheet. 


(A3) 
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The velocity W . H is computed by multiplying the input values W/W 

mxa, ni • j ___ 

? mia, m 

(44) and W cr m (27). With equation (A3) the midchannel velocities at the hub and tip 
can then be calculated from the value of W mid m - For computing the integration in 
equation (A3) , the interval r^ to r^ is divided into eight equal intervals . Since the 
integration is started at the mean radius, only four integration steps are required in each 
direction. The integration is done numerically using the trapezoidal rule. In the case 
where /3 is given as input only at hub, mean, and tip (rather than at all nine stations) 
linear interpolation is used to determine the values which are not given. The meridional 


streamline curvature l/r„ is assumed to vary linearly from hub to tip. 

/ W W 

The ratios of surface velocity to midchannel velocity | — and — E- 


can be 


W. 


'mid mid/ 

computed based on an assumption of linear variation of either streamline curvature or 
streamline radius of curvature. These equations are derived in the next section. If 
linear variation of curvature is assumed, 


n r 


W 


= e 


\8 nj 8 1 


1 

2 


£,c . 


(A4) 


W 


mid 


where n is the distance along the orthogonal from the blade suction surface. Thus, the 
velocity can be determined at any point on the orthogonal from (A4). If, on the other 
hand, linear variation of radius of curvature is assumed. 


W 


W 


•= < 


mid 



< c s - V S- 


n 

n oJ 


c + c 

P S 


n 0 C s C p 

V C s 


-4 -a 


“ V c s 


(A5) 


^ C =C e 
P s 


The blade surface velocity at hub, mean, and tip can be computed by substituting n = 0 
and n = n Q in equation (A4) or equation (A5) and using the results from equation (A3). 
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With the velocities obtained, the weight flow past the orthogonal can be computed 
from 

»calc ‘ill fo° PW in dr <A6> 

The inner integral is computed at hub, mean, and tip using trapezoidal integration over 
eight equal intervals from suction to pressure surface. 

If JY = 1 or 2 (linear curvature or radius of curvature variation), the critical veloc- 
ity ratios W/W are computed using the input value of W and equations (A4) or (A 5), 
and (A3). Finally, the weight flow parameter pW/p M W„„ is calculated from 


pW 

p”W 
H cr 


1 - 2 ^-! 
y + 1 



(A7) 


If JY = 3 (linear variation of static pressure), the blade surface static pressures 
are calculated based on the surface velocities from 


JL = 


1 - 


y - 1 


w 


y 

nr- 1 


r + l \w 


cr , 


(A 8) 


and pressures along the orthogonal are calculated from 


p _ p s + n / p p p s 


p M p" 


n. 


0\ P 


,, pM 


(A9) 


Finally 



1 


- 

r-i~ 


pw _ 

@1 

y + 1 

1 _ 



pW" cr 1 

y - 1 

A 1 

w 







J 


1/2 


(A10) 


Note that the blade surface velocities are based on linear variation of either curvature 
(eq. (A4)) or radius of curvature (eq. (A5)) and the velocity used to calculate 
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pW/p”W (eq. (A10)) is based on linear variation of static pressure. This is not con- 
sistent. 

After computing the inner integral in equation (A6) at hub, mean, and tip, the outer 
integral is approximated by Simpson’s rule. If JZ = 1, no further calculations are made. 


If JZ = 2, further calculations are made with new values of 
the weight flow w as a function of ^W/W cr j mid 


( W /W cr ) 


mid, m 


to get 


The value of 


( W /*cr) 


m 


giving maximum weight flow is estimated using para 


bolic approximations 
than 0. 001, the computation is stopped 


mid, m 

When two successive estimates of 


(*/*«) 


mid, m 


differ by less 


The maximum weight flow is normally deter- 


mined within five iterations. 

If JZ = 3 or 4, further calculations are made to determine a value of ^W/W cr ^ 


mid, m 

- . If w . is less 
giv giv 

which will give 


that will give a weight flow w which is close to the input value w ( 
than choking weight flow, there are two values of ^W/W cr ) mi( j m 
w calc = w giv‘ ^ JZ = 3, the smaller, or subsonic, value of m will be 

found. If JZ = 4, the larger, or supersonic, solution will be 'found. When 


| w . ■ w calcl 001 w giv’ ca l cu l a ti° ns are stopped. A solution is normally found 

within five iterations. If w^ iy is larger than the choking weight flow no solution exists. 
In this case calculations are made for five values of (w/W cr ^ m j d m , and the choking 
weight flow can usually be estimated from these values. 


Derivation of B lade -to-B lade Velocity Variation (eq. (A5)) 

The method of calculating blade surface velocities from midchannel velocities is 
based on reference 2. The assumptions that the flow is steady relative to the blade, 
nonviscous, and isentropic along the blade to blade streamline orthogonal yield 


dW _ W 
dn r c 


(All) 


This is equation (1) of reference 2 and can be derived from the force equation. 

If it is assumed that the curvature (C = l/r c ) varies linearly along the orthogonal, 
the equation can be integrated to obtain 


n r 


W 


2(C - C ) 

' r'v o ' 


W 


= e 


mid 



(A12) 
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which is equation (7) of reference 2. Since the curvature is assumed to vary linearly, 

c = c s + (C p - C s ) -2- (A 13) 

n 0 

When equation (A13) is substituted in equation (A12), equation (A4) is obtained. 

If it is assumed that the radius of curvature varies linearly, then 

r c = (r c>s + [< r c >p-< r c>s]^ 


(A14) 


Using this in equation (All) we have 


dW ; 
dn 


W 


< r c ) s + [ (r c'p‘ <r c ) s]f- 


(A15) 


Integrating equation (A15) from the midchannel gives 


log 


W 


n„ 


W 


mid 


^ r c^p “ ^ r c^s 


log; 


r 

< r c>s + 

blp. 

- (r c>s 


2 

(r c>s + [ ( 

r c>p - 

(r c»Jr 

0. 




(A 16) 


When C g / C p , equation (A16) can be solved for W/W mid after substituting r c = 1/C, 
yielding 


n«C C 


W 


W mid 


2 l C p + ( C s * C p> 


C p +C s 



(A17) 


For the special case C = C , equation (All) may be integrated (r is constant), and 


s p 


the resulting equation solved for W/W m ^ d , yielding 


W 


n °\2 " n of s 


W mid 


= e 


(A 18) 


This completes the derivation of equation (A5). 
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cttd F|v INPUT DATA SHEET 

The lowest line in each block displays Roman numerals for card sequence and arabic numerals for card 
columns within each card. 

[All numerical values (columns ® - ©, ©) are entered in nine-column fields, (read in with a 
FORTRAN format specification of E9. 5) usually as ±xxxxx±zz representing ±.xxxxx multiplied by 


1 


I 

V40 (MAY BE ALPHANUMERIC) 


0 


II 

1 

n 2-40 (MAY BE ALPHANUMERIC) 




q 


q 

(5) 




( 14 ) 

IMI 

“0” 

[MI 

lM_ 

JX 

JY 

JZ 

KR4 

KR5 

KR6 

n 

KR7 

r 

l/sec 

u 

lb/ sec 
w giv 

i/rt 

(l/r C ) h 

i/ft 

(UrcH 

ft 

r h 

ft 

r t 

/ 

/ j 

z 

/ 

/ j 

J* 1 


/.V 

ytf. 

o.$<*n 

O 

<3 

t>j7f 

/ 

III 

1 

1 

2 

— i 

4 

5 

6 

7 

IV 

1-9 

10-18 

19-27 

V 

1-9 

10-18 

19-27 

28-36 


O’) 

~(20) 


C 22 ) 

"MI 


(25) 

~W~ 

MI 

[MZ 

ft 

1/ft 

1/ft 

ft/ sec 

lb/ft 2 sec 

ft 

i/ft 

1/ft 

ft/ sec 

Itfft 2 sec 

n O,h 

c s,h 

c p, h 

w cr,h 

P"Vh 

n O, m 

c s, m 

c p, m 

W 

¥V cr, m 

P"Vm 

o.cdxof 

/$" 7/ 


■7V7. 7 


t.cs+5 

,4/t> 

g¥-Z 

937-4 

4-9*53. 

VI-A 

1-9 

10-18 

19-27 

28-36 

37-45 

VI-B 

1-9 

10-18 

19-27 

28-36 

37-45 




(31) 


M 

ft 

n 0,t 

i/ft 

c s,t 

1/ft 

s.t 

ft/ sec 
w cr,t 

ll/ft 2 sec 
P"W C r,t 

07 tS 

tZ 60 

6.30 

C> 

A? . A" V 

VI -c 
1-9 

10-18 

19-27 

28-36 

37-45 



1*T 

"W 



O 9 ) 

isr 

C 41 ) 

(42) 

deg 









-5:^ 

-/6. 3 








VII -A 
1-9 

10-18 

19-27 

VII-B 

1-9 

10-18 

19-27 

VII-C 

1-9 

10-18 

19-27 



(43) 



z station 
code 

^cr^mid, m 

0 

fooo 0003 

o.t 

VIII 

1 

2-9 must be 
(NUMERIC) 

10-18 


Figure 5. - input data sheet. 
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SAMPLE PROBLEM 
NEGATIVE REACTION TURBINE 
( 12 ) GAMMA = O. 14000E 01@ OMEGA 


0.48800E 03 


HUB 

MEAN 

TIP 




N ZERO 


C SUB S 


= 0. 59I70E 00 

W SUB CR 


5 0.420806-01 £ 

3 0- 1571 OE 

02 

^ 0.882006 

01 ( 
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-0.10650E 02 
-0.19350E 02 
^ -0. 28500E 02 


VSH* 0.13422E 01 VSM= 0.141456 01 VST* 0.14594E 01 VPH* 0.801066 00 VPM* 0.76449E 00 VPT* 0. 75868E 00 


IK 

LITTLE A 

LITTLE 8 


ALPHA 

1 

-0.868136-02 

-0.85064E 

02 

-0.20347E-03 

2 

-0. 201076-01 

-0.132886 

03 

-0. 47 l 26E-03 

3 

-0. 352566-01 

-0. 18037E 

03 

-0. 82632E-03 

4 

-0.534636-01 

-0.227436 

03 

-0.125306-02 

5 

-0.741406-01 

-0.27393E 

03 

-0. 173776-02 

6 

-0.989616-01 

-0.32339E 

03 

-0.231946-02 

7 

-0.125596 00 

-0.371926 

03 

-0. 29435E-02 

8 

-0.153496 00 

-0.41941E 

03 

-0. 359736-02 

9 

-0. 182 146 00 

-0.46571E 

03 

-0. 42690E-02 

IK 

M( IK) 

N( IK) 



1 

0. 704246-02 

-0.844676 

02 


2 

0.63676E- 02 

-0. 13204E 

03 


3 

0. 50701 E-02 

-0. 17946E 

03 


4 

0. 299076-02 

-0.226756 

03 


5 

0 . 

-0.273936 

03 


6 

-0.405716-02 

-0.32470E 

03 


7 

-0. 93200E-02 

-0.375416 

03 


8 

-0. 158616-01 

-0.42612E 

03 


9 

-0.237276-01 

-0.476896 

03 



PH= 0.33632E 02 PM* 0. 


PT = -0* 70389E 02 OH* 0.10071E 01 QM = 0.10000E 01 QT = 0.976556 00 


WSH * 0.10422E 01 HMH = 0.7T652E 00 WPH * 0.62204E 00 
WSM = 0. 11316E 01 WMM * 0.80000E 00 WPM = 0.61159E 00 
WST = 0. 12331E 01 WMT * 0.844926 00 WPT = 0.64103E 00 
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00 
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00 

0.624636 

00 

0.633896 

00 

0.62903E 

00 

0.61557E 

00 

0.597406 

00 

0.577176 

00 

0.556616 

00 

0.53685E 

00 


0.80000E 00 WT FLOW CALC 


0.61444E 00 


-0.384396-01 


WSH = 0.97410E 00 WMH * 0. 72576E 00 WPH = 0.581386 00 

WSM = 0.10609E 01 wNH = 0. 75000E 00 WPM = 0.57337E 00 

WST = 0.11625E 01 WMT = 0.79656E 00 WPT = 0.804346 00 


A( IK ) = 

0.63343E 

00 

0.6261 3E 

00 

0.61270E 

00 

0.595666 

00 

0.576816 

OO 

0.557396 

00 

0.53826E 

00 

0.519996 

00 

0.502936 

00 

A ( IK ) = 

0.631146 

00 

0.63304E 

00 

0.623546 

00 

0.506956 

oo 

0.586386 

00 

0.563996 

00 

0.541226 

00 

0.51904E 

00 

0. 49802E 

00 

A ( IK ) = 

0.61421E 

00 

0.63222E 

00 

0.632046 

00 

0.620466 

00 

0.60237E 

00 

0.581106 

00 

O.55805E 

00 

0.537036 

00 

0.516536 

00 


0.750006 00 WT FLOW CALC 


0 .601276 00 0 = -0.16173E-01 


WSH = 0.905976 00 WMH * 0.675006 00 
WSM = 0.990136 00 WMM * 0.700006 00 
WST = 0. 10919E 01 WMT = 0.748196 00 


WPH = 0.54072E 00 

WPM = 0. 5 351 46 00 

WPT = 0.567646 00 


A( IK ) = 

0.627196 

00 

0.61340E 

00 

0.595356 

00 

0. 5 7509 E 

00 

0.554066 

00 

0.53323E 

00 

0.51324E 

00 

0.494506 

00 

0.477236 

00 

A ( IK » = 

0.63386E 

00 

0.626496 

00 

0.610406 

00 

0.589296 

00 

0.565726 

00 

0.541426 

00 

0.517546 

00 

0.49477E 

00 

0.473556 

00 

A ( I K ) = 
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00 

0.633746 

00 

0.625106 

00 
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00 
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00 

0.561926 

00 

0.53812E 

00 

0.515436 

00 

0.494476 

00 


0.700006 00 WT FLOW CALC = 0.584016 00 0 = 0.129906-01 


WSH = 0.93436E 00 WMH 

WSM = 0.101966 01 WMM 

WST = 0. 112146 01 WMT 


0.69615E 00 
0. 72083E 00 
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WPH = 0.55766E 00 

WPM = 0. 551 0 76 00 

WPT = 0.582936 00 
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00 

0.61943E 

00 
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00 

0. 5841 8E 

00 

O.56390E 

00 

0.54367E 

00 
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00 

0.505406 

00 
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00 

A{ IK ) = 

0.6 336 5E 

00 

0.62999E 

00 

0.61651E 

00 

0.59717E 

00 

0.57476E 

00 

0.55U9E 

00 

0.527716 

00 

0.505146 

00 

0.48397E 

00 

A ( IK )* 

0.622876 

00 

0.63385E 

00 

0.628616 

00 

0.61352E 

00 

0.593 06 E 

00 
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00 
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00 
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00 
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00 
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Figure 6. - Output listing. 
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